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ABSTRACT 

We relate the star formajtion from cold baryons condensing in virialized structures 
to the X-ray properties of the associated diffuse, hot baryonic component. Our compu- 
tations use the standard "semi-analytic" models to include and connect three sectors 
of the complex astrophysics involved: first, the formation of dark matter halos through 
accretion and merging, after the standard hierarchical clustering; second, the star for- 
mation governed, after the current "recipes", by radiative cooling and by feedback of 
the supernova energy into the hot baryonic component; third, and novel, the hydro- and 
thermodynamics of the hot phase, rendered with our Punctuated Equilibria model. So 
we relate the X-ray observables concerning the intra-cluster medium (namely, the lumi- 
nosity - temperature relation, the luminosity functions, the source counts) to the thermal 
energy of the gas pre-heated and expelled by supernovae following star formation, and 
then accreted during the subsequent merging events. 

Our main results are as follows. At fluxes fainter than Fx ~ 10^^^ erg/cm^ s the 
X-ray counts of extended extragalactic sources (as well as the faint end of the luminosity 
function, their contribution to the soft X-ray background, and the Lx — T correlation 
at the group scales) increase considerably when the star formation rate is enhanced for 
z > 1 as indicated by growing optical/infrared evidence. Specifically, the counts in the 
range 0.5-2 keV are increased by factors ~ 4 when the the feedback is decreased and star 
formation is enhanced as to yield a fiat shape of the star formation rate for 2 < z < A. 

Such faint fluxes are well within the reach of next generation X-ray observatories 
like AXAF and XMM. So very faint X-ray counts will soon constitute a new means of 
gaining information about the stellar processes (formation, and supernova feedback) at 
z > 2, and a new way to advance the understanding of the galaxy formation. 
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1. INTRODUCTION 

The history and the fate of cosmic baryons con- 
stitTitc an OTitstanding problem in cosmology. It is 
matter of current debate their time-shifting partition 
among the stars and the diffuse cold and hot compo- 
nents, the latter being heated by stellar activity or at 
the expenses of the gravitational energy in the dark 
matter (DM) condensations. 

Two lines of evidence enter the argument: one 
comes from the canonical bands, optical-UV and IR; 
the other will come - we argue - from the X-ray band. 

The first assessments from 0-UV data of the global 
star formation rate (SFR) suggested a peak at ^ « 1.5 
with a sharp decline out to 2: w 4 (Madau et al. 1996; 
Connolly et al. 1998); in such a picture, less than 20 % 
of the stars would have been formed at z > 2. Based 
on the hierarchical cold dark matter (CDM) models 
of structure formation, which predict the galaxy as- 
sembly to be a gradual process proceeding through 
merging and accretion events, such history has been 
modeled (sec White & Frcnk 1991; Cole et al. 1994; 
Kauffmann, White & Guiderdoni 1993; Baugh et al. 
1998) by strongly coupling the star formation (SF) 
inside the halos to their dynamical growth; in fact, 
the SFRs have been assumed to scale up considerably 
with the DM halo masses. 

However, such scenario has been seriously chal- 
lenged by more recent observations (Pettini et al. 
1997; Dickinson 1998; Meurer et al. 1997) indicat- 
ing that the apparent decline for z > 2 was seriously 
affected by dust extinction; the underlying SFR may 
be increased by factors 3-15, up to yield a flat plateau 
for 2; > 2 (Madau, Pozzetti & Dickinson 1998). This 
picture is confirmed by the large statistics obtained 
in denser fields by Steidel et al. (1998). Support- 
ing indipendent evidence comes from the still scanty 
IR data obtained with SCUBA (see Small, Ivison & 
Blain 1997; Hughes et al. 1998; Barger et al. 1999). 

But the issue is still at stake. In fact, a flat SFR 
might easily lead to an overproduction of the IR back- 
ground, and of the metal adundance in high-z ab- 
sorbers (see Ellis 1998). Moreover, it is not clear how 
the optical fields sampled so far should be weighted 
to yield a reliable estimate of the average SFR. 

Luckily, there is another angle to the baryon his- 
tory. This is provided by the hot baryonic phase 
which is related to the SF and can be observed di- 
rectly in X-rays as intra-cluster medium (ICM). A 
large amount of baryons corresponding to fift/O w 0.1 



is observed at z « in galaxy groups and clusters on 
scales around i? ~ 1 Mpc (White et al. 1993; White 
& Fabian 1995; see also Fukugita, Hogan fc Peebles 
1998). Such baryons are heated up to virial temper- 
atures kT ~ GM/IQR in the gravitational potential 
wells corresponding to masses M w 10^^ — 10^^ M©; 
at densities n « lO^'' cm"'' they emit by ther- 
mal bremsstrahlung copious X-rays with luminosities 
LocR^v? TV2 « 1044 erg/s. But the ICM is not se- 
cluded; indeed, the shape of the Lx — T correlation, 
which is observed to be bent from L oc in very 
rich clusters (Allen & Fabian 1998) to a much steeper 
shape in groups (Ponman et al. 1996), indicates that 
the temperature of the external, infalling gas plays a 
key role in determining the X-ray properties of small 
groups and poor clusters (Cavaliere, Menci & Tozzi 
1997, 1998; Balogh, Babul & Patton 1998). Such ex- 
ternal temperature, in turn, may be set by supernova 
feedback to values exceeding the virial value in small 
merging lumps, so a connection is envisaged between 
the stellar processes and the astrophysics of the X-ray 
emitting ICM. 

This connection will provide a relationship be- 
tween the O, UV or IR observations and the X-ray 
data, which is the scope of this paper. In Cava- 
liere, Menci & Tozzi (1997, 1998a, 1998b hereafter 
CMT97, CMT98a, CMT98b) we have already devel- 
oped a semi-analytical approach to the X-ray emis- 
sion, to deal with the DM merging histories, with 
the associated infall of external preheated gas, and 
with the resulting shock compression, all concurring 
to estabilish a sequence of approximate hydrostatic 
equilibria of the ICM (Punctuated Equilibria model). 
Our predicted, bent Lx — T correlation is in eccellent 
agreement with the observations; the parameter of 
the model is the initial temperature of external gas, 
which was assumed at a constant value around 0.5 
keV adopted from simple estimates in the literature. 

Such a temperature is of key importance as it hap- 
pens to be close to the virial temperatures in typ- 
ical groups, and may be surmised to rise gradually 
during structure formation. In addition, the amount 
of gas expelled by supernova feedback is essential in 
determining the amount of baryons filling of shallow 
potential wells, and hence the X-ray luminosity of 
small groups. So here we set out to compute these 
processes in full, based on the detailed recipes cur- 
rently implemented in the semi-analytical models for 
star and galaxy formation (SAM hereafter, see Kauff- 
mann, White & Guiderdoni 1993; Cole et al. 1994; 
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Baugh ct al. 1998; Somerville & Primack 1998) in 
hierachically evolving DM halos. This allows us to 
connect in a single framework the hierarchical clus- 
tering, the astrophysics of star forming baryons, and 
the thermo- and hydrodynamics of the X-ray emitting 
plasma. In particular, we will show how forthcoming 
observations with AXAF can X-ray the SF history, so 
complementing the 0-IR information. 

2. THE MODEL 

Toward our scope wc need to relate both the SF 
and the astrophysics of the ICM to the gravitationally 
dominant DM content of the halos. After the hierar- 
chical clustering paradigm, the latter form due to the 
gravitational instability of initial random density fluc- 
tuations, and then evolve through stochastic merging 
of smaller or sometimes comparable units into larger 
structures. The resulting statistics is well estabilished 
in the form of the extended Press & Schechter theory 
(EPST, sec Bond ct al. 1991; White & Prenk 1991; 
Bower 1991; Laccy & Cole 1993). 

As to the process of SF, we adopt an analytical ren- 
dition of the SAM approach. In the above papers the 
baryons gravitationally bound to a DM halo, after be- 
ing heated to its virial temperature, are assumed to be 
exchanged among three phases: the cold phase pro- 
duced by radiative cooling; the resulting condensed 
phase into stars; the hot phase constituted by gas 
preheated by supernova explosions and further raised 
to the virial temperature of the potential wells. The 
model includes coalescence of baryons following halo 
merging, and the luminosity-color evolution due to 
the rise and fall of successive star generations. 

We include the X-ray emission in the very same 
framework, using the model proposed by CMT97: at 
each time step the X-ray luminosity is computed from 
the hot phase of the ICM, and from the boundary 
conditions relating this to the external gas infalling 
in each merging event. As shown in CMT97-98, such 
boundary conditions are essential in providing a real- 
istic model for the X-ray emission from clusters and 
groups. In fact, a set of equilibrium states exists 
for the ICM in a given halo, each corresponding to 
a different merging history and hence to a different 
boundary condition. The average X-ray luminosity 
corresponding to a given mass M is then obtained 
by convolving over all possible merging histories the 
emission corresponding to each equilibrium state. 

Because of the added complexity to the standard 



sense, and of the related increase in computing time, 
we choose to express the merging histories analyti- 
cally with the EPST rather than using the equivalent 
Montec Carlo simulations. Below we describe the ba- 
sic ingredients of our model. 

Whenever necessary, and untill otherwise specified, 
the redshift z will be associated with the epoch t on 
adopting a critical Universe with Ho = 50 km/s Mpc. 

2.1. Gas Cooling, Feedback and Star Forma- 
tion Rates 

Our recipes to relate the SFR to the dynamics of 
DM halos are basically taken from Cole et al. (1994) 
and are summarized as follows. 

a) For each DM mass M (or circular velocity Vc oc 
M^/"^), the baryonic content is divided into a cold gas 
phase [apt to radiatively cool within the halo] with 
mass nic, into a star phase with mass m*, and into 
the complementary hot phase with mass m/j. The 
initial gas content associated to all the halos sums 
up to the universal baryonic density fif,. The initial 
stellar content at Zi ~ 10 is taken to be nil. 

b) The mass of the cold phase is increased from 
inside out by cooling processes: at each time step At, 
we have then Atoc(i'c) = '^t^t^qqi PgArcooi, where gas 
density Pg and the dependencies of the cooling radius 
rcooi on Vc and on t are given, e.g., by Somerville & 
Primack (1999). We adopt the cooling function used 
by Cole et al. (1994), which was computed assum- 
ing a primordial mixture of 77 % hydrogen and 23 % 
helium. 

c) Then it is computed the amount of cold gas 
which goes into stars, and the corresponding mass 
is transferred from the cold phase to the star phase. 
Such amount Ami,{vc, t) = m* At is regulated by the 

SFR, namely rh,*(uc,i) = mc/Tsr. The timescale is 
parametrized in the form t^{vc) = t° (i;c/300 km/s)"* ; 
the cold mass mc is updated after each time step by 
subtracting not only the mass of stars which have 
formed but also the cold gas expelled given below. 

d) A final transfer is due to supernova feedback 
from the cold back to the hot phase; this involves the 
amount Amt = /3{vc) Am*(wc,t), where the feedback 
per unit star mass is given by /3 = {vc/vh)^°''^ , with 
the parameters Vh and ah- 

For each circular velocity Vc, we compiite the fur- 
ther variations of the masses of the gas components 
and of the star content due to merging of DM ha- 
los (treated in detail in §2.3), possibly followed by 
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galaxy coalc;sccncc. The latter is included on con- 
sidering the probability that the galaxy coalescence 
time (parametrized as in Cole et al. 1994) is shorter 
than the halo survival time (depending on Vc), and 
by averaging over all possible halo merging histories 
(details are given in Poll et al. 1999). The values of 
the parameters r^, a*, Vh and ah are discussed in §3. 

2.2. The X-ray Emitting ICM 

At the time t, the X-ray luminosity due to bremsstrahlung 
emission by the hot baryons m/j in a halo of mass M 
is given by (CMT98b) 



Af^) G\M,M')T^{M)p'/^I{M,t). (1) 



Here r^(M) = 4{M/1Q^^ Mq^/^ (1 + z) keV is the 
virial temperature corresponding to M; the average 
DM density inside the cluster is p oc (1 + z)^; the nor- 
malization constant A will be adjusted so as to match 
the height of the observed local Lx — T relation at 
T = 4 keV; the shape factor / describes the internal 
ICM distribution. Finally, G is the density jump (ra- 
tio) across the shock induced at the cluster boundary, 
at around the virial radius, by the infalling gas; this 
will depend not only on the cluster mass M, but also 
on the mass M' of the infalling clumps. 

The contribution of emission lines has been in- 
cluded using the standard and public Raymond-Smith 
code. 

We stress why the Lx — T dependence differs from 
the self-similar power-law L cc T^. To a small degree 
this is due to the shape factor I{M, t), a slowly vary- 
ing function of M and t of the detailed form given by 
CMT98b. This includes the integration over the clus- 
ter volume of the internal density run (normalized to 
the value at the virial radius); the latter is provided by 
the hydrostatic equilibrium once the gravitational po- 
tential (/>(M) associated to M is given (we shall adopt 
the form given by Navarro, Frenk & White 1997). 

Much more important is the strong dependence 
of on the halo merging histories; in fact, the 
(squared) compression ratio G^ at the shock connects 
the ICM with the infalling gas associated with the 
merging partner of mass M'. Of this, the fraction /* 
hovers around the clump M' having being expelled 
and heated to a temperature T^{M') by supernova 
feedback; the complementary fraction 1 — /* is still 
contained inside the virial radius of M' at the virial 
temperature Ty{M'). During a merging event, both 



components fall into our well of mass Af, and the 
compounded G^ entering the luminosity (1) reads 

G\M,M') = Ug^[T,XM')] + {l-U)g^[T,{M')]. (2) 

Here g{T) denotes the shock compression factor for 
the inclusion of each gas fraction into a cluster or 
group, given in CMT98a,b to read 



5(T) = 2(l- J) + [4(l 



T \ 2 
T2. 



T 

T2 



1/2 



(3) 



the expression for T2 is given by the same authors, 
and for strong shocks reads simply 

kT2{M) « ^mH<t){M)/3 + 3kT/2 , 

where tuh is the proton mass. As anticipated above 
eq. (2), the argument T in eq. (3) takes one of the 
following two values: either the virial temperature 
(with the coefficient given by Metzler & Evrard 1997) 

n {M') = 4 (M'/lO^^ M0)2/3 (1 + z) kcV (4) 

for the 1 — /* baryons retained in the lump M' ; or the 
temperature of the ejected baryons provided (in the 
timestep Af) by supernova feedback. 



kT.{M')= EsNVSNTh, At 



(5) 



where Esn = 10^^ erg/s is the energy per super- 
nova, r}sN is the number of supernovae per solar mass 
(3.210-3 for the Scalo IMF we shall use here); the 
latter equality holds after point d in Sect. 2.1. Note 
that the ratio m*At/Am/j = {vh/vc)°"' depends only 
on the current circular velocity of the halo and not on 
the progenitor masses. 

Finally, the cold gas fraction /, expelled outside 
the virial radius of M' and heated at T* (M') is com- 
puted adopting the Cole et al. (1994) recipe, i.e., 
assuming that the reheated gas is expelled from the 
halo. So the fraction of gas reheated in the timestep 
At reads 

A/* = Anih/mh ■ (6) 

Note that the quantity A/* depends on the merging 
history of the halo, as will be described in detail in 
Sect. 2.3. 

Our guideline will be that, since the X-ray luminos- 
ity in eq. (1) depends strongly on /, and which in 
turn depend on the star formation and on feedback, a 
strong correlation must exist between the X-ray emis- 
sion and the SFR. 
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2.3. Convolutions over the Merging Histo- 
ries 

Both the SF and the X-ray emission depend ulti- 
mately on the merging liistory of the halos. In fact, 
at each merging civent the gas reservoir of a halo are 
increased by the gas (cold and hot) enclosed in the 
merging partner, and this changes the processes of 
cooling, feedback and SF together; in addition, by the 
sequence of overlapping merging events, the infalling 
gas compresses the ICM at the cluster boundary, and 
sets the X-ray luminosity (eq. 1) through the factor 

So to obtain the observable quantities - average 
value and scatter - associated to a halo, we must con- 
volve over all merging histories leading to that mass. 
In detail, we adopt the following procedure, essen- 
tially an analitical rendition of the SAM, but comple- 
mented with the proper description of the physics of 
the X-ray emitting ICM. 

i) We define a time-velocity grid, with grid size 
Ai = io/100 and Avc = 10 km/s {Iq is the present 
cosmic time). 

The number density N{M, Zi) of DM halos with 
mass M [corresponding to a circular velocity Vc = 
(10Gff(z)M)i/3] at the initial redshift ~ 10 is 
taken from the Press & Schechter (1974) expression. 
Initially, to each halo we associate a galaxy with the 
same Vc of the halo (more than one galaxy per halo 
is an exceedingly rare circumstance for Zi ~ 10). For 
each Vc the corresponding mass M and virial tem- 
perature Ty are computed. As said above, the initial 
gas content associated to all halos corresponds to the 
universal baryonic density fi;,, and the initial stellar 
content is taken to be nil. 

For each circular velocity Vc, the baryonic content 
is divided into stars, cold gas and hot gas as described 
in §2.1 under point a). 

ii) At the next time step, we compute the mass 
transfers: hot cold stars — > hot, as described in 
§2.1, point b), c) and d). 

iii) For each circular velocity Vc, we compute the 
further variations of the mass of the gas components 
and of the stars content due to merging of DM ha- 
los, possibly followed by galaxy coalescence. We also 
compute the compression factor (eqs. 2) and the 
corresponding X-ray luminosity (eq. 1). 

Since the merging process is stochastic, at each 
time step and for each Vc (corresponding to M) we 



compute the probabihty density d'^P{M', t\M,)/dM'dt 
that a halo of mass M' has merged with a halo of com- 
plementary mass M — M,' to yield the considered M 
in the time interval At; within the standard hierar- 
chical clustering, such a probability is provided by the 
EPST. 

The average mass of the cold and the hot gas con- 
tents in a halo with circular velocity Vc are updated 
according to following equation, analogous to White 
& Prenk's (1991): 



mh{vc,t + At) 



ruhivct) + At J^'-^'^Um' X 

N(M',t) d^P{M',t\M) 

N{M,t) dM'dt "W^^ ) ■ \l ) 



This yields the average increment due to the merging 
together with the halos of previous baryon reservoirs. 
An analogous equation holds for the baryonic mass 
in stars m.,{vc,t) and in the cold phase mc{vc,t), as 
well as for the fraction /* of ejected gas, whose time 
increment is given by eq. (6). 

As for the X-ray luminosity, this is computed from 
the compression factor G^(M, M') (eq. 2) inserted 
into eq. (1). Its average value is computed by an 
analogous convolution to the above one: 



Lx{vc,t + At)=Lx{vc,t) + At dM'x 

■h) 



fM{vc) 

1 

N(M\t) (fl-'{M'.I^M) 
N(M,t) dM'dL ^ 

A (^i^yc^M, M') T^{M)p'/^{z){8) 



iv) The values of rric, m* and m/, are reset for every 
Vc after all their (positive or negative) increments due 
to cooling, SF, supernova feedback and merging have 
been computed in the steps ii) and iii). 

v) Finally, for each Vc the associated optical lumi- 
nosity at the wavelength A is computed by the convo- 
lution 



Sx 



Jo 



t')m^{t')dt' 



(9) 



where the spectral energy distribution of luminous 
stars 4'\{t) is taken from a canonical model of stel- 
lar population synthesis in its latest version (Bruzual 
& Chariot 1998). 

vi) We increment the current time by At, and re- 
peat the whole procedure through steps ii) to v) until 
the output time is reached. 

We comment that our analytical rendition of the 
galaxy formation sector is technically similar to the 
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formulation by White & Frcnk (1991), but differs in 
a number of respects: a) we consider here aU three 
components (stars, cold gas phase, hot gas phase) 
which arc involved in the SF process, while in the 
paper by White & Frenk (1991) only cold gas and 
stars where considered; b) we implement the compre- 
hensive recipe used in SAM for the star formation and 
feedback, while in Frenk & White (1991) a simplify- 
ing assumption of self-regulation had been adopted to 
obtain a simple expression for the SFR. To compute 
the optical luminosities we have used here a standard 
population synthesis model in the updated version by 
Bruzual & Chariot (1998). On the other hand, we 
plan to add chemical enrichment models in an exten- 
sion of this work. 

Our main improvement concerns of course the in- 
sertion in the old framework of the novel issues con- 
cerning on the X-ray emitting baryons. 

3. CHOICE OF PARAMETERS 

In fact, our aim is to obtain the X-ray counterpart 
of different SF histories. We consider two extreme 
cases A and B, both consistent with the present data 
discussed in the Introduction. 

Model A is characterized by a SFR declining be- 
yond 2; « 2, as was originally suggested by Madau et 
al. (1996), and as was obtained from the SAM (Baugh 
et al. 1998). From the latter authors, we adopt the 
set of parameters reproduced in the first row of our 
Table 1, and chosen by them to yield an acceptable fit 
to the local galaxy luminosity function. The resulting 
SFR we find (see fig. la) is peaked at 2; » 1.5 and in 
fact declines considerably for z w 2. 

The corresponding galaxy luminosity function at 
2; = is shown in fig. lb. As discussed by Cole et al. 
(1994, see their Table 1; also Baugh et al. 1998), both 
the decline of the SFR at large z and the flat shape of 
the galaxy luminosity function are due to the strong 
supernovae feedback at small Vc resulting from their 
set of star formation and feedback parameters. The 
same feedback aff'ects strongly the X-ray luminosities 
through the compression factor and the factor m\ 
in eq. 1. 

Model B instead is characterized by a SFR flat be- 
yond 2; ~ 2, see fig. 2a. The set of parameters leading 
to this are shown in the second row of Table 1. Such 
choice corresponds to a milder feedback even in small 
halos, consistent with recent works by Ferrara & Tol- 
stoy (1999) and by Martin et al. (1999). This yields a 



rather steep local galaxy luminosity function (fig. 2b) 
below L«, a feature considered tenable, if marginally, 
in view of the data by Zucca et al. (1997). 

In deriving the galaxy luminosity functions we fol- 
lowed Cole et al. (1994) in normalizing the mass to 
light ration in terms of T, the total mass in stars 
divided by the mass in luminous stars with mass 
> 0.1 Mq. We adopted T = 2.7 for Model A (the 
value of the fiducial model by Cole et al. 1994), and 
T = 2 for Model B. 

For comparison with the SAMs, we adopt thir 
"fiducial" cosmological/cosmogonical context, namely 
the CDM model with ag = 0.67, O = 1, Ob = 0.06 
and h = 0.5. In all cases, the Scalo IMF has been 
adopted, as widely used in SAMs. 

4. RESULTS 

The outcomes in X-rays from Model A are shown in 
fig. 3. The Lx — T relation provides a good fit to the 
data {panel a) down to the sub-keV range, the region 
most affected by supernova feedback (eq. 6). Note 
the mild evolution out to z — 1. The corresponding 
X-ray luminosity function {panel b) also agrees well 
with the local data (Ebeling et al. 1997; De Grandi 
et al. 1998) and shows little evolution out to z « 1 
in agreement with recent data (Rosati et al. 1998). 
We also plot {panel c) the complementary contribu- 
tion to the soft X-ray background (details are given 
in CMT98a), which is well below the observational 
limits even when a 75 % resolved contribution is sub- 
tracted out of the ROSAT data (see Hasinger et al. 
1998; Giommi, Fiore & Perri 1998). 

The corresponding results for the opposite extreme 
Model B are shown in fig. 4. Here the Lx—T relation 
is flatter at low T, and constitutes a marginal flt to 
the existing data. This is because the lower amount 
of gas expelled by virtue of the smaller feedback (with 
a minor effect from the increased temperature T*) 
reduces the changes with T of the compression fac- 
tor G entering the X-ray luminosity (eq. 1); the re- 
sult is closer to the gravitational self-similar scaling 
L (xT^. This goes back to the circumstance (see fig. 
2) that, in the absence of an appreciable contribution 
/, (7^(T*) from supernovae in the merging partner M', 
the prevailing term (1 — /*) 5^(71,) is determined only 
by the virial temperature ratio Ty{M)/Ty{M'). Such 
ratio of purely gravitational temperatures does not 
change appreciably with T when averaged with the 
self-similar merging rates; so a nearly constant 
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obtains, and eq. (1) tends to L oc T^. 

Another feature of Model B is the increase with z 
of the X-ray luminosities; this is due to the larger 
amount of hot gas [the factor (to^/M^) in eq. 1] 
made available by the stronger global, supernova ac- 
tivity, but retained in shallow wells by virtue of the 
weak feedback. This has a number of implications: 
not only the faster evolution of the Lx — T relation 
represented by the dashed curve in panel a, but also 
the increased number of faint sources represented in 
the luminosity functions N(L,z) in panel b; in addi- 
tion, the contribution to the soft X-ray background is 
larger, and barely consistent with the observational 
limits. 

But the key X-ray test telling apart the two models 
is provided by the source counts A''(> F) shown in fig. 
5. The larger number of sources predicted at high z 
by model B implies faint counts larger by factors ~ 
4 relative to Model A; their redshift distributions are 
shown in fig. 6 A and 6B. Correspondingly, in fig. 7 we 
show for the two models the luminosity density in X- 
rays, the counterpart for the hot baryons of the SFR 
associated with condensed baryons shown in fig. 1. 
Note that for model B the peak is higher and shifted 
to ^ w 2, corresponding to the larger SF activity at 
high z. We also stress that, in any case, the effect 
of SF on the X-ray emission is always delayed by the 
few dynamical times taken by the merging activity 
(see cq. 8) to affect the ICM. 

Finally, we show in fig. 8 the correlation between 
the X-ray flux and the optical {B and K) magnitudes 
in DM halos, as predicted by Model A and by Model 
B. 

5. CONCLUSIONS 

In this paper we have related two quantities: the 
X-ray emission from the hot diffuse baryons in hi- 
erarchically evolving potential wells, and the history 
of the baryons condensed into stars. We have pre- 
sented a computational package which grafts on the 
extended Press & Schechter treatement of the hierar- 
chical clustering both our model for the thermo- and 
the hydrodynamics of the X-ray emitting ICM, and 
our analytical implementation of the current recipes 
adopted in the semi-analytical models of galaxy for- 
mation. 

We have shown that supernova feedback is an es- 
sential ingredient in deriving not only the early SFR, 
but also a realistic shape of the Lx — T correlation 



for groups and clusters of galaxies, and in predict- 
ing counts and luminosity functions of faint X-ray 
sources. In fact, the SFR at 2; ^ 2 depends strongly 
on the amount and the temperature of the baryons 
retained is small potential wells. But as such small 
halos are included into larger structures at lower z, 
the same quantities affect the X-ray properties of the 
accreting halos. 

Specifically, we find that if the SFR is peaked at 
z w 1.5 with a deeline to higher z, then the local 
Lx—T correlation flattens strongly going from groups 
to rich clusters of galaxies, and evolves little with 
z; the corresponding luminosity density in X-rays is 
peaked at 2; ~ 1. Conversely, if the SFR was already 
high since 2: « 4, then a smoother local Lx — T rela- 
tion obtains, closer to the self-similar form Lx oc 
down to poor groups; relatedly, the counts of such 
sources brighter than Fx — 10~^^ erg/cm^ s exceeed 
the former case by a factor ss 4 in the energy band 

0.5. Even fainter fluxes are within the reach of the 
next generation X-ray observatories like AXAF (and 
XMM, though near its confusion limit). To be conser- 
vative with surface brightness, we have included in the 
counts only Lx > erg/s. In addition, the above 
excess counts has been computed with a very con- 
servative low energy threshold (0.5 keV) for AXAF, 
which in fact may have an effective threshold of 0.25 
keV (R. Giacconi, private communication), while 0.1 
keV is planned for XMM; for example, on lowering 
the threshold down to 0.25 keV the excess would be 
doubled (see fig. 5). 

All that will make in the near future the X-ray ob- 
servations of groups and clusters a powerful means to 
discriminate the strength of the feedback processes, 
which determine also the SF history at high z. For 
example, the excess X-ray counts by a factor 4 cor- 
responds to a filling by stellar and hot baryons of 
the early potential wells which yields a SFR larger by 
~ 10 at z w 3 - 4. 

We have investigated the sensitivity of our results 
to changes of the cosmological framework. We find 
that the difference between Model A and B persists 
when one considers different cosmologies within the 
usual constraints provided, e.g., by the cosmic age and 
the local counts of X-ray clusters. In particular, we 
have checked that for a flat A-cosmology {flo = 0.3, 
= 0.7, fib = 0.04) the count excess is still larger 
than 3.5. A more extended study of the SFR- X-ray 
connection for different cosmologies will be presented 
elsewhere. 
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We have explored the effect of varying the star for- 
mation and feedback parameters, and of changing the 
star formation recipe altogether. We have first kept 
a* fixed at 0, and tried values of intermediate be- 
tween Model A and B, yielding intermediate decline of 
the SFR at 2; ^ 2. We find X-ray counts again exceed- 
ing those of Model A for any value of 2 < < 5.5; 
e.g., at Fx = 10~^^ erg/cm^ s for = 3.5 the excess 
over Model A is a factor 2.3, while for ah = ^ the 
excess is still factor 1.9; recall that when the energy 
band is extended down to 0.25 keV the difference be- 
tween this case and Model A is doubled. Thus, X-ray 
observations with low energy threshold, as planned 
for the next space observatories, can pinpoint the ef- 
fective strength of the feedback within the range de- 
limited by our extreme Models A and B. If we now 
change a» to the other extreme value -1.5 the above 
results are not changed appreciably. This is because 
it is the feedback which dominates (see Cole ct al. 
1994) the amount and the thermal state of all baryons 
in shallow, early potential wells, and so governs both 
the X-ray emission and the SFR at large-z. 

We have also investigated the effect of adopting 
recipes for star formation and for feedback different 
from Cole et al. (1994). In particular, we imple- 
mented in detail the recipe by Kauffman et al. (1998), 
where all the gas is retained in the halos and the feed- 
back parameter has about the same Uc-dependence of 
our Model B. This yields a flat SFR at ^ > 2; but 
in addition we obtain X-ray counts close to those in 
Model B, confirming the trend: more stars - more X- 
rays. However, the lower normalization of the feed- 
back leads to a very flat local Lx — T correlation. 

Finally, we have checked that our results are not 
sensitive to the inclusion of additional heating or cool- 
ing. Specifically, the supernova heating of the hot 
fraction inside the deep wells turns out to be not rel- 
evant, unless the associated T, is forced \vp to values 
close to the virial temperature around 10 kcV, at vari- 
ance with the accepted star and supernova energetics. 
Such findings are consistent with those by Mctzler 
& Evrard (1997) based on hydrodynamical/N-body 
simulations. On the other hand, considering explic- 
itly the cooling of the gas ejected beyond the virial R, 
the temperature decreases by only 5 — 10%, due 
to the low densities of the gas in the external regions. 
This holds both for a distribution spread outside R 
with a density following the DM's, and for a constant 
density shell of thickness v^sc located at R. 

In sum, the thermal state and amount of hot 



baryons (with their X-ray emission) and of the con- 
densed baryons in shallow potential wells (yielding the 
SFR at high z observed in the 0-UV bands) both de- 
pend on the same energy feedback from supernovae. 
The X-ray and the 0-UV bands concur to yield pre- 
dictions and information concerning the baryon his- 
tory; in particular, the X-rays catch directly the feed- 
back in the act and probe its effective strength. We 
add that, while in this paper we have focussed on 
the truly hot component (at T ^ 10^ — 10^ K, mak- 
ing up a minor but vocal fraction of all baryons), the 
connection here investigated involves other baryonic 
phases, observable in different bands; e.g., the luke- 
warm (T < 5 10^ K) gas contributing to the the dif- 
ferent SF histories can be compared with the amount 
of photoionized gas probed by absorbtion in the Ly-a 
clouds (both in and out the galaxies). Such a variety 
of independent probes will stricly bound the amount 
and behaviour of the feedback processes, which are 
presently highly uncertain but in fact govern the early 
cosmic SFR and so the emission properties of the faint 
galaxy population. 
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FIGURE CAPTIONS 



Fig. 1. — Panel a): The global star formation rate, 
as a function of redshift z for Model A. 
Panel b): The B-band luminosity function of galaxies 
in Model A; for comparison, we show the data by 
Zucca et al. (1998). 

Fig. 2. — Same as fig. 1, but form Model B. 

Fig. 3. — Results from Model A: 
Panel a): The Lx —T correlation at 2 = (solid line) 
and at z=l (dashed line). Group data from Ponman 
et al. (1996, solid squares); cluster data from Marke- 

vitch (1998, open triangles). 

Panel b): The local (solid line) and the z = \ (dashed 
line) luminosity function. For comparison we show 
the data by Ebchng et al. (1997). 
Panel c): The contribution of hot baryons to the 
soft X-ray background for sources with Fx < 410~^^ 
erg/cm^s, is compared with the total observed values 
(open stars, Hasinger et al. 1997), and with the 25% 
unresolved component (see Giommi, Fiore, & Perri 
1998, solid squares). 

Fig. 4. — Same as fig. 2, but for Model B 

Fig. 5. — The source counts from Model A (solid 
line) and from Model B (dashed line), in the energy 
band 0.5-2 keV. The dashed region corresponds to the 
observed ROSAT cluster number counts by Rosati et 
al. (1998). 

Fig. 6. — The redshift distribution of X-ray sources 
is plotted for differents fluxes: 10~^^ < F < 10"^'' 
erg/cm^ s (solid line), and the very faint range 
10"^^ < F < IQ^^'' crg/cm^ s (dashed line). Panel 
A: Model A; panel B: Model B. 

Fig. 7. — The X-ray luminosity density (for photon 
energies in the range 0.1-10 keV) as a function of red- 
shift. Panel A: model A; panel B: model B. These 
constitute the X-ray counterparts of the 0-UV lumi- 
nosity density corresponding to the SFR shown in fig. 
1 

Fig. 8. — The correlation between X-ray fluxes and 
optical magnitudes is plotted (shaded region) for B 
magnitudes (upper two panels) and if-magnitudes 
(lower two panels), for Model A (left column) and 
Model B (right column) 
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Table 1 
Model Parameters 





' * 


a* 


Vh 




Model A 


2.8 Gyr 


-1.5 


140 km/s 


5.5 


Model B 


2.8 Gyr 





140 km/s 


1.5 



The parameters corresponding to the two reference SFRs introduced in the 
text. These define, as a function of the halo circular velocity Vc, the star formation 
time scale r*(t>c) = r° (fc/300 km/s)"* and the mass Aruh = {yc/vh)~°''^ Am*(t'c,t) 
reheated by SN feedback. 
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I 

O 

a 




1 2 3 4 5 6 

z 




-23 -22 -21 -20 -19 -18 -17 -16 -15 -14 -13 

Mg-5 1og(h) 








1 2 3 4 5 6 

z 




-23 -22 -21 -20 -19 -18 -17 -16 -15 -14 -13 

Mg-5 1og(h) 




-2-1012 

log (Lx/lC^erg s-') 
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-2-1012 

log (Lx/lC^erg s-') 





Fy [0.5 — 2 keV] erg cm ^ s-1 



